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We consider a broad class of interactions between radiation and a light scalar field, including both 
conformal and disformal couplings. Such a scalar field potentially acts on cosmological scales as 
dark energy and could also appear in modified gravity theories. We study the consequences of these 
couplings on the mixing between the scalar field and photons in galaxy clusters in the presence of 
a magnetic field. In particular we focus on the resulting turbulence-induced irregularities in the 
X-ray and UV bands. We find new bounds on the photon-to-scalar couplings, both conformal and 
disformal, which complement laboratory experiments and other astrophysical constraints. 
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I. INTRODUCTION 

Light scalar fields may play a crucial role in the recent history of the Universe. They could be the seeds for the 
acceleration of the expansion of the Universe and play the role of dark energy at late time [T] . They could modify the 
properties of gravity on very large scales [5]. Their observable properties depend crucially on the type of couplings 
they have with matter fields and radiation. Conformal couplings involve a Yukawa-like interaction between scalars 
and fermions. Their study is therefore very motivated from a fundamental perspective. They are responsible for 
the modification of gravity by scalar fields and explicitly imply the changes of matter trajectories which are tightly 
constrained by solar system experiments |3]. On large cosmological scales, the effects of a modification of gravity by 
conformal couplings could be seen in the growth of structure and will be actively sought for by the Euclid mission 
[1]. These conformal couplings do not lead to a direct interaction of scalars to photons at tree level as the photon 
Lagrangian is conformally invariant classically. Such a coupling is only induced at the one loop level when virtual 
fermions mediate the interactions between two photons and one scalar . We will include such a conformal coupling 
in our analysis. Another type of coupling is particularly important and will be one of the original aspect of this article: 
disformal couplings mil]. They involve the coupling of scalars to the full energy momentum tensor of matter and 
radiation. They could arise if gravity becomes massive and the scalar field arises as the zero helicity part of a massive 
graviton [S] . These couplings involve two derivatives of the scalar field and can therefore only play a role in dynamical 
situations. They have no effect on static tests of gravity involving a modification of static forces for instance and are 
therefore not testable in the solar system (Sj- On the other hand, such a coupling to photons has several consequences. 
They are very tightly constrained by astrophysics [10] , collider experiments and cosmology [111 [12] . Astrophysically, 
the coupling to baryons would increase the burning rate of stars and supernovae. Similarly, the primakoff process 
in stars induced by the disformal coupling to photons implies that a photon in the electric field of a nucleus would 
release two scalars. This effect allows one to obtain constraints on the suppression scale of the disformal coupling to 
photons m- In collider experiments, two quarks would transform into two scalars which would appear as missing 
energy. Both LHC experiments ATLAS and CMS were able to constraint the coupling scale too m- The disformal 
coupling plays also a role on large scales where it would induce a change of the speed of light which would distort 
the cosmic microwave background spectrum m- In the laboratory, in the presence of a magnetic field, the photons 
of a laser beam would see their polarisation rotate and become elliptic in a way akin to the well-known effects for 
Axion-Like-Particles (ALP). In this case, the disformal coupling together with the direct coupling to photons lead to 
a non-zero transition probability between an initial photon state and a scalar [T^. The disformal coupling alters the 
usual transition probability for ALP’s. In this article, we will concentrate on the effects of both the conformal and the 
disformal couplings in astrophysical situations where the external magnetic field appears inside large galaxy clusters 
for instance. 

In galaxy clusters, the magnetic fields are turbulent, as revealed by radio observations. The idea here is to observe 
a source of radiation (hereafter X-ray photons and UV light) that lies at the center of a galaxy cluster. In the regime 
where the photons and the scalar field mix at the source, the content of the radiation along the line of sight is actually 
a mix of photons and scalar field excitations. Because of the turbulent nature of the magnetic field and the energy- 
dependent mixing, strong irregularities can be introduced in the photon energy spectra, in a way that would resemble 
a set of random absorption lines. The exact shape of the induced spectral shape is unpredictable, but the statistical 
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properties of the induced noise can be inferred m- This is used to constraint parameters of the scalar field model in 
a similar way to the ALP parameter space (mass and coupling to photons) which has been analysed using gamma-ray 
sources m, and X-ray sources [T7]. 

In this paper, we consider the electrodynamics of disformally coupled scalars in section 11. In section III, we apply 
it to the propagation through turbulent media. In section IV, we deduce constraints on the parameter space of scalars 
coupled to photons from X-ray observations, and similarly in section V with UV probes. In section VI, we obtain 
constraints on the parameter space of the scalar models. We conclude in section VII. 


II. DISFORMALLY COUPLED SCALAR FIELDS 
A. The scalar models 

We consider the coupling of a light scalar field to matter governed by the action 

s = J ( l67rGjv ~ ~ + Smii/Ji, Qf,^) + , (1) 

where Gn is Newton’s constant, which is related to the reduced Planck scale as mpj = The first term in 

the action is the Einstein-Hilbert action describing General Relativity (GR). The scalar field is canonically normalised 
with an interaction potential V{(j)) which we do not specify here. In our study we will take the potential term to 
be a simple mass term V{(f>) = This simplified setting could serve as a template for massive gravity models 

where the scalar polarisation of the massive graviton becomes a low mass scalar field. More complex cases could be 
considered in models of screened modified gravity where the potential has a minimum (j) which could depend on the 
environment, i.e. it would depend on the matter density. Expanding around this minimum cj) = cj) + S(j), the potential 
would reduce, to leading order, to the case of a massive particle corresponding to the field 6(j). In the physical situations 
that we will consider where photons traverse large clusters in the Universe, we will assume that the scalar field in this 
environment reduces to a massive scalar with a quadratic potential. The matter action Sm describes the matter part 
of the standard model of particles physics. The coupling to photons in directly involves the scalar field. Indeed, 
this coupling between the scalar and photons is dictated by the disformal metric 
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~ T ^4 dfi(f)du(j ), (2) 

which will constrained using galaxy cluster observations in the following sections. We have not considered a disformal 
coupling to matter particles in Sm as it must be heavily suppressed to comply with the ATLAS constraint. Indeed it 
is shown in [101 that the dimensionful disformal coupling to baryons is constrained at the level of Mb > 490 GeV. As 
a result, in this paper, we consider that at the astrophysical energies (X-rays and UV) considered here the disformal 
coupling to matter is negligible and will be set to zero. Notice that the coupling to different matter species together 
with the one to photons does not have to be universal. Moreover, the astrophysical and collider constraints obtained 
in | 10 | indicate that disformal couplings have very different effects in different contexts, from laboratory experiments 
to the burning of stars, and that they are better analysed as separate constants whose values should be inferred from 
either experiments or observations. In this paper, we focus on the disformal coupling to photons. 

The metric is the Jordan frame metric of photons with respect to which photons are conserved = 0 , 

where the Jordan frame energy momentum tensor is . On the other hand, the metric 5 ^ 1 / defines the 

Einstein frame and in this frame energy-momentum is not conserved. The Klein-Gordon equation is modified by the 
presence of the energy momentum of photons 

= ( 3 ) 

The non-conservation of photons, which results from the coupling of the scalar to photons and therefore the exchange 
of energy-momentum between the scalar and radiation, appears in the second term of this equation. In what follows 
we will restrict ourselves to the leading order effects of the disformal coupling between the scalar field and photons. 
Therefore we calculate only to leading order in 1/M^, and to this order we have where is the 

energy-momentum of radiation in the Einstein frame. Higher order terms in 1/M^ are present but they would lead 
to correction terms in 1/M® in the Klein-Gordon equation that we shall neglect. For radiation, this implies that we 
will use the energy momentum tensor 

_ y_ Tp2 

(7) “ “ 4 


( 4 ) 




3 


The Klein-Gordon equation reduces to 


= ( 5 ) 

where we will explicitly take V{(j)) = ^rri^cjP. At this order, the action can be expanded as 

S = j ( iGt^Gat ~ -y{4>) + . (6) 

These models have links with varying speed of light theories which can be found in m- The coupling scale M 
is constant and unknown and should be fixed by observations as we will see. In certain models where the screening 
effects appear due to non-linearities in the kinetic terms, the action for the canonically normalised fluctuations 54i 
around the background value (j) is such that the coupling scale M becomes environment dependent. One may therefore 
consider that M may depend on the local density. For this reason, we shall quote our new bounds by recalling that 
they are obtained from galaxy clusters which are much less dense than stars. As a consequence, the star burning 
constraints may be relaxed, for instance, in models of the K-mouflage type where stars are screened whereas clusters 
are not, implying that the suppression scale in stars and in clusters can be very different. A thorough study of the 
density dependence of the scale M is left for future work. Here we only consider the specific environment provided 
by galaxy clusters from which sound bounds can be extracted. 

In the following, we will also specialise to the case where the new scale M = yjmmp\. This is motivated by the 
decoupling limit of the new theories of massive gravity where the scalar field appears as the scalar polarisation of a 
massive graviton of mass m. In this case the disformal coupling is parameterised by the scale M related to the mass 
of the graviton m. The burning of stars such as the sun can be increased by the disformal coupling to photons as 
it can lead to a new Primakoff effect whereby scalar escape stars and affect the life-time of stars. This leads to the 
astrophysical bound M > 346 MeV, valid in dense environment around 100 g/cm^. As already noted, this bound 
may be relaxed in models where the disformal coupling scale is density dependent. If this is not the case and when 
= mmp\, this implies that one cannot probe masses lower than 10“^® GeV which corresponds to scales of order 
less than 1 km. We will find that the constraints obtained from galaxy clusters where the density is way lower, i.e. 
10“^“^ g/cm^, are looser and masses as low as 10“^^ GeV can be allowed. 


B. Electrodynamics with a disformal coupling 

Let us generalise the previous setting by introducing a field dependent coupling constant, controlled by a new 
unknown scale A, so that the kinetic term for photons contains both the disformal coupling and a direct coupling 

'S'rad ^ - J (^1 + > (7) 

where contractions are made with the Jordan frame metric. The fine structure constant becomes field dependent 
«{(/)) = —^ where a is its value in the absence of the direct coupling. The purely disformal case is obtained by 

taking A —)■ oo. When the scalar field has a time dependence at the background level 0(t), this leads to a time 
variation of the fine structure constant a(0) which is constrained at the 10“® level at small redshift z < 3 [18]. In the 
following we will consider that the background evolution (j) of the scalar field (p is such that this constraint is satisfied. 
Denoting by 0o the value of the background field now, we identify the fine structure now as a = a{(j)o). The variation 
of the fine structure constant is given by 


a A J- 4(/) 

which is small provided that \p — pol A in the recent past of the Universe. In the following, we will use the 
simplifying assumption that p A. Our results can be easily extended to the case when this is not the case anymore. 

As mentioned before, the disformal coupling leads to a coupling between the scalar field and photons which becomes, 
to leading order in 1/M^, 




(9) 
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This is the action that we shall use to study the mixing between the scalar and photons. The equation of motion for 
the photon gives the generalised form of Maxwell’s eauation|14| 


dn, 


4(j) 


1 




^a/3 


- - d^<I^F^)] = 0 ■ 


( 10 ) 


where the leading term in XjM^ has been kept. In the following we consider that there is a background magnetic field 
Bi and that we decompose the vector field into 


Ai — c^jkFjXfi^ T ( 


( 11 ) 


and Aq = 0. We work in the Lorentz gauge dfj^A^ = 0 which implies that, in Fourier modes, the vector field ai is 
transverse kia^ = 0 with respect to the wave vector ki. Similarly we assume that the scalar field takes a background 
value ([> that we shift 


(j) = (l) + 5(j) 


( 12 ) 


where the fluctuations S(f> are massive with a mass m. We also assume that the time variation of the background field 
(j) is negligible in the astrophysical situations we consider, e.g. the propagation of photons through a galaxy cluster. 
Defining the dimensionless constants 




we find that Maxwell’s equation becomes 


2 _ L _ _ 

“ “ A ’ ^ “ M4 


4i 


(1 - a ){E - k )ai - -^CijkBjkkSfj) = 0 


(13) 


(14) 


where and E is the energy of the photon represented by the wave ai ^ k.x)^ Klein-Gordon 

equation becomes 


(1 + - k^) - + 2b^ k^ - kik 


B^B^ 

' B2 


4^ 

dcj) -\- jCi}^ — 0 


(15) 


where all the contracted indices involve the Kronecker 5ij. It is convenient to specialise the propagation along the z 
axis and normalise 


54> —>■ vT+If ai —>■ \/l + a? ai. 


(16) 


The mixing between the photons and the scalar field is such that the dispersion relation is still to leading order E ^ k, 
and upon using k = idz we have 


{E - idz)V + MV = 0 


(17) 


where we have defined the vector 


and the mixing matrix 


M = 


O; 

V = \ a. 


y 

iS(j) 


0 

0 


0 

0 

2Bx 


__ 

AV(l+62)(l-a2) 

2Bx 

Tie Cl _ 


V ~ AyM+bAA-a'^) Av'(l+f'D(l-aD ~ 2B(T+b2) “T I+52 
In a plasma the mixing matrix is modified and becomes 

/ “ol 2B. 


M = 


"'pi 

2E 

0 

2B„ 


"'pi 

2E 


AA/(l+f)2)(l-a2) 

2B, 


Av^(l+f) 2 )(l-a 2 ) 


2B^ 


b^E 


V ^^/(l+^Aj{VAA) A^y(l+b^){l-a^) 2B(l+b2) A- 


B2 - 


(l-|t) 


(18) 


(19) 


( 20 ) 
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where uj^i = is the plasma frequency. Notice that the plasma frequency is not field dependent as we have 

assumed that 0/A ^ 1. This allows one to deduce the the probability of conversion from a photon to the scalar in a 
uniform magnetic field after a distance s. It can be exactly calculated and reads 


with 



1 




sin^ 


27rs 


Api 

Aa 

Ab 

A 


2E 

b^Eil-§) 

2E{l + b‘^) 1 + 62 

2B 

aVtTF 

47r 

V(Apl - Aa)2 


( 21 ) 


( 22 ) 

(23) 

(24) 

(25) 


When 6 —0, the ALP case is retrieved. In the following we will see that the presence of the new coupling 6 introduces 
drastic differences in the energy dependence of the transition probability. 


C. Energy dependence 


We are interested in cases where the transition from photons to scalars is not suppressed. This happens when the 
prefactor — (a |-a close to unity. The conversion is therefore efficient when we have that Ab > |Api — Aa|/2. 


1 + 


As Api and Aa depend on the energy E of the photon, only a band of energies satisfies this condition. More precisely, 
there is one term proportional to the energy. 


and one inversely proportional 


A^ = 


b'^E 


A I, = - 


2E 2E(l + b^)' 


(26) 


(27) 


We have neglected the component and put a <C 1 to comply with the bounds on the variation of a in the recent 
past of the Universe. These terms are represented on Fig. [^for three different configurations. We have denoted by 
the energy at which A^/^ crosses Ab. When E^ < E^, an efficient conversion takes place in this energy range. 
If not, no efficient conversion occurs. In the middle panel of Fig. [^the case when E^, ~ E^ (resonant conversion) is 
also illustrated. 

We have that 

^ ^ |m2-a.g^(l + 6^)|A 

8bVTT¥ 

^ ABy/TT¥ 

~ A 62 ■ 


(28) 

(29) 


For b = 0, El is equal to the critical energy in the ALP case above which the mixing between photons and ALP’s is 
efficient. The A^ and A^ terms cross at Ex = i/iAyE/ . We can now find the range of parameters m and 6 for which 
an efficient conversion occurs, that is to say Ei <C E^. For this, we solve Ei = E^ for 6 and get 


6 


2 

X 


Hi 


± Vd 


(30) 
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Efficient conversion between | and Resonant conversion No conversion 



Ei log(E) Ipg(E) Ipgjgj E^ 


FIG. 1: Illustration of the different conversion regimes depending on the critical energies: E_i < E^: efficient conversion 
between E^ and E-f (left panel), E^ ^ Ef-. resonant conversion (middle panel), E^ > E-f-. no efficient conversion, the interval 
of energy where conversion would be possible does not exist anymore (right panel). 


where we have defined 


D = 




128- 


4i^ 


A 2 


(31) 


Two cases need to be distinguished. When < ujpi{l + &^), we have that bx is unique and equal to 


b 


2 

0 


32^2 


(32) 


and an efficient conversion is only possible when b < bo- Here the conversion is only efficient in an interval in energy 
between and and this is only possible for small values of the ratio b = B/M^ < bo. This is a new case compared 
to the usual ALP situation as the conversion there only involves one threshold energy and it becomes efficient above 
this threshold. When + 6^), we have two solutions for bx given by 


b 


2 

1 


32^2 
A^m^ ’ 



(33) 


In this case, when b < bi the conversion is efficient while for larger values of b it is only efficient in the narrow band 

bl<b^ < ^ -1. 

^pi 

In some model realisations, it is possible to relate b and m using either the ansatz b = B/mA or for massive gravity 
models b — B/mmpi. We have summarised the results using these relations in Fig. where we have assumed values 
for B and of 10 and 0.01 cm“^ respectively that are typical of galaxy clusters, and a conformal coupling scales 
A = 10^ GeV (left) and A = lO^i GeV (right). In Fig. § the blue plain curves separate the two domains, where above 
this limit, no efficient conversion occurs. The blue dashed line is the region where resonant coupling happen, when 
6 ^ psj rfi— 1, in which case —>■ 0. This requires a big fine-tuning of b since the width is very narrow. The red lines 

n 

are just two examples of specific cases for b. For both panels of Fig. the red plain curve relates b and m using the 
scale A and the red dashed curve using wpi. Finally we have represented the horizontal like 6 iini = B/Mn^ ~ 10“^® 
where Miim = 346 MeV is the limiting value allowed by the burning rate of stars. Only values of b < &iini are allowed 
when the constraint on M is density independent. When M and A are related using = toA, it appears that only 
a small mass window is accessible, just below the m = ojpi line. If that constraint is relaxed, the mass range where 
signal is present lies somewhere between ~ 10“^® eV and 10“^^ eV for A = 10^ GeV, and between ~ 10“^° eV and 
10“^^ eV for A = 10^^ GeV. So considering larger conformal scales should narrow down the accessible mass window. 

Note that when the relation b = B/mmp\ is considered, the comparison between our potential limits and the 
ones from the star burning rate should be taken with caution, the scale M should depend on the density because 
of screening effects. When the suppression scale M depends on the density, such as for K-moufiage models, the 
limiting behaviour provided by the stellar burning constraint can be relaxed in a model dependent way. Indeed the 
density inside clusters is typically p ^ 10 “^^ g cm“^ whereas it is typically p ^ 100 g cm“^ in stars and even a mild 
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A = 10^ GeV 



m (eV) 


■ Coupling limit 


- Resonant coupling 


h = 


h = 


B 

mA 

B 

mmpi 



m (eV) 

- Coupling limit - b = 

- Resonant coupling - 


B 

mA 

B 

mmpi 


FIG. 2: Different regimes in the parameter space b versus m for values of A = 10^ GeV and A = 10^^ GeV (right), see text 
for details. 


dependence on the density may render the star constraint weaker than the bound obtained from galaxy clusters. The 
analysis of this dependence is left for future work. 

In the following, two energy bands will be considered for the photons: X rays between 1 keV and 7 keV, and UV 
around 1 eV. From the above analysis of the energy dependence of the signal, it is already possible to estimate the 
ranges in A and b that will be accessible to this study. This is illustrated in Fig. § where A is fixed at 10^ GeV (left 
panel) and 10^^ GeV (right panel). The region where an irregularity signal is potentially observable {E^ < is 
highlighted in green. The two energy ranges for photon observation are displayed as horizontal zones. From the left 
panel of Fig. |^we infer that (if high enough) the signal will be accessible for b lower than ~ 10“^^ in UV and ^ 10“^^ 
in X rays. Increasing A up to 10^^ GeV narrows the signal region and make it inaccessible to both energy ranges. It 
is therefore expected that it will not be possible to reach higher A values. 

In a similar fashion, one can fix values of b to draw accessible regions in the (U, A) plane. This is illustrated in 
Fig. where the signal regions are drawn for two values of b, 10 (left) and 10 (right). From these figures, it 
appears that in the case when b — 10“^^, accessible values for A are below 10^ GeV and 10^ GeV for X rays and UV 
respectively. When b is reduced to 10“^^, these values shift to 10^^ and 10^ respectively. 


D. Laboratory experiments 


The disformal coupling to radiation has effects on the the propagation of laser beams in optical cavities in the 
presence of a transverse magnetic field. Indeed, the mixing between the photons and the scalar field implies that the 
polarisation of an initially linearly polarised beam rotates and becomes slightly elliptical. The strongest constraints 
are given by the PVLAS experiment m which cannot test models when A > 10^ GeV. Better constraints are obtained 
from the light shining through a wall experiment at DESY m where the probability than an initial laser beam is 
regenerated on the other side of a wall due to the crossing of the wall by the scalar field. In this case, the experiment 
cannot constrain models with A > 10^ GeV. For lower values of A, the masses of the scalar field are extremely 
constrained when = mpim, with m < Hq ~ 10“^^ eV M- All these laboratory tests have been performed in 
densities p ^ 10“^^ g.cm“^. This should be compared with the available mass range obtained using b < bo 
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FIG. 3: Illustration of the accessible values for b, 

6 = 10 “^^ 

— E-^ 



FIG. 4: Illustration of the accessible values 


for A = 10’’ GeV (left) and A = 10“ GeV (right). 


b = 10"^^ 



and = m mp\ leading to 
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(34) 
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which can be as low as 10 GeV. There is also the astrophysical limit b < &ii„i = 

M,? 

m > (35) 

mpi 

The latter prevents to probe masses of order Hq and in fact m > 10“^® GeV. Hence, laboratory constraints imply that 
astrophysical tests of the disformal coupling can yield new result only for large values of A > 10^ GeV and masses 
m > sup(m*,rniim). This can be seen in Figj^where m\\^ > rrii, for A = 10^^ GeV. In summary, the limits that we 
will derive from galaxy clusters cannot probe the very small mass range where m ~ Hq but larger masses which are 
still very small compared to the masses of known particles. 


III. PROPAGATION IN TURBULENT MAGNETIC FIELDS 

The case of a homogeneous magnetic field that was considered in the previous section is a simplified picture that 
is not representative of astrophysical magnetic fields. For instance, measurements of the Faraday rotation of the 
polarization of synchrotron emission from radiogalaxies embedded in galaxy clusters showed that magnetic fields in 
clusters are usually turbulent on various spatial scales m, following a Kolmogorov power spectrum {P{k) oc k 
up to scales of tens of kpc |22lE^ . The Milky Way also hosts a turbulent magnetic field, in addition to the large 
scale ordered component (see |25| for a recent model), that is found to a have a Kolmogorov power spectrum on 
scales from 0.3 pc to 100 pc pS]. Notice that the interaction between the scalar field and photons do not affect these 
measurements, which are in the radio band. 

In the following, the turbulence of the magnetic field is described as in m by summing over a large number 
of plane waves having randomly chosen propagation directions. The plane waves have different wavenumbers and 
corresponding amplitudes so as to describe the power spectrum assumed for the turbulence. This method has the 
advantage of generating a three-dimensional turbulent magnetic field that can be computed at any point in space 
coordinates, based only on a random set of angles for the propagation direction and phase of each plane wave. While 
one realization of the turbulent magnetic field corresponds to a random choice of every angles, it is obvious that the 
exact 3D realization of the turbulent magnetic field cannot be measured in nature and thus remains unknown. In 
the description of the turbulence with a Kolmogorov power spectrum, the fluctuations on the smallest spatial scales 
are less important than the fluctuations at the highest scale and therefore become irrelevant as the scales decrease. 
The line of sight for the propagation in the turbulent magnetic field of the system, initially consisting of photons, is 
then divided in small domains such that in each domain the magnetic field is assumed to be homogeneous, such that 
equation [T7| is easily integrated. 

The initial beam is conservatively assumed to be unpolarized. This is accounted for in the density matrix formalism: 

p = \ Ay \ ® Ay icj))* , (36) 

and an initial density matrix po = diag(Y 2 , V 27 O). After propagation in every magnetic domains, the density matrix 
is computed recursively as 


Pk+i = e 


— iMk s 


• Pk 




(37) 


where s is the size of the domain crossed by the beam and Aikis the mixing matrix given by Eq. |20[ The 3-dimensional 
magnetic field used in the mixing matrix is computed from the turbulence model randomly generated and is assumed 
to be homogeneous in each domain. The value and direction of the magnetic field between the domains are related 
through the power spectrum assumed in the turbulence model. 

It has been shown in m for the ALP case, that the photon survival probability after crossing various magnetic 
domains, that can be computed with Eq. |37| has a very complex energy dependence around the critical energy. This 
effect is due to the turbulent nature of the magnetic field. The irregular pattern in the photon survival probability 
around the critical energy depends on the exact realization of its turbulence. In the present case where the mixing 
only occurs for energies in the energy range between and spectral irregularities are present around both 
critical energies. This signature is similar to the ALP case with the addition of a new parameter b, that modifies 
the irregularities pattern when b ^ b^. For this reason, constraints that have already been obtained on ALPs with 
the spectral irregularities signature (see usiini) cannot be translated directly to constrain disformal couplings. As 
explained in Sec. we are mainly interested in constraining the very low mass domain for the scalar field. This can 
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be achieved with the addition of the plasma term in the mixi ng m atrix, as it was considered in |17| . so that the 
constraints are valid for all masses m ^ ujpi ~ 10“^^ eV (see Eq. 281. 

An example of the photon survival probability from an initial unpolarized beam of pure photons is shown on Fig[^ 
for three values of b and A assuming propagation in a rando m rea lization of a magnetic field typical of galaxy clusters 
with B = 1/rG and Ue = 0.01 cm“^. As shown in section IIC 
range for the strong mixing is reduced. 


when b is increased, decreases and the energy 
Nevertheless, the expression of shows that this could be in principal 


compensated by lower values of A. For lower values of A, the oscillation length of the mixing is decreased and the 
spectral irregularities are tighter. This is shown on the second panel of Fig. where A = 10^° GeV as compared to 
the first panel where A = 10^^ GeV. For this reason, two methods are proposed in the following. The first method 
has been presented in [12] and is based on X-rays of the Hydra A galaxy cluster. Due to the limited energy resolution 
of the instrument, of the order of 5%, this method can probe spectral irregularities at high A, and corresponding low 
b. This case is outlined on panel 1 and 2 of Fig. where the gray area shows the sensitive energy range. In order to 
identify spectral irregularities on small energy scales, a spectral resolution of the order of 10“^ is required. A second 
method based on high resolution spectroscopy of a BL Lac object, PKS 2155-304, is then used to obtain constraints 
on lower values of A but higher values of b (see panel 3 of Fig. |^. 



10'^ 10'^ 10"' 1 10 10^ lO'^ 

Energy [ keV ] 


FIG. 5: Photon survival probability for three values of 6, assuming propagation in one realization of the turbulent magnetic 
field model of Sec |IV| The grey area shows the energy range sensitive to X-ray observations and the grey line the energy domain 
probed by the UV spectroscopy of PKS 2155-304. 


IV. CONSTRAINTS FROM X-RAY OBSERVATIONS 

Because of typical values for the electron density and magnetic fields in galaxy clusters, the critical energy E^ is 
expected to be of order of a keV. Gonsequently, spectral irregularities due to the conformal mixing are expected to 
occur in this energy band, as shown in Fig. Spectral irregularities should then be searched in the spectrum of a 
bright X-ray source embedded in a galaxy clusters with a strong magnetic field. As explained in HZ], the non-thermal 
emission from the central AGN of Hydra A is a good target for the present search because of the strong magnetic 
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field in its cluster that has been extensively studied in the literature 


A. Magnetic fields in the Hydra A galaxy cluster 

The magnetic field in the Hydra A galaxy cluster can be probed with Faraday rotation maps from the radio emission 
of the lobes of the central AGN. To obtain a measurement of the magnetic field, a profile of the electron density is 
assumed. It has been estimated in pS] by using the X-ray surface brightness measured by ROSAT in m with the 
method of |32| . The electron density profile is parameterized as a function of the radius r from the central source 
following : 


Heir) = 



(38) 


with nei = 0.056 cm“^, ne 2 = 0.0063 cm“^, Vci = 33.3 kpc, rc 2 = 169 kpc and /3 = 0.766. The same electron density 
profile is used later in the simulations of the propagation of the photon/scalar system to compute the plasma term 
along the line of sight. Moreover, the r.m.s. of the turbulent magnetic field at a radius r in the cluster is believed to 
scale with the electron density as B{r) oc ne(r)“^ |33| . 

The measurement of the magnetic field also depends on the geometry that is assumed for the orientation of the jets 
and lobes from which the Faraday rotation maps are produced. The angle of projection of the northern jet on the 
line of sight 0 explains the depolarization asymmetry that is observed between the two lobes [2H]. The uncertainty on 
this parameter still remains large, as 0 is constrained in the range between 30° and 60°. In the latest analysis [23], a 
magnetic field strength Bq at r = 0 of 21 ^G is found if 9 = 30°, compared to Bq = 36 /xG if 6* = 45° and Bq = 85 /xG 
if 9 = 60°. To be conservative, Bq = 21/xG is assumed in the following. In the same analysis, a most likely value 
of 1 for the scale parameter is found. This value is also used in the present work. For the turbulence of the 
magnetic field, a Kolmogorov power spectrum is measured on spatial scales up to 10 kpc. For the simulations of 
the turbulent magnetic field performed in this study, a minimal scale of 1 kpc is assumed. As explained in im, 
spectral irregularities rapidly become irrelevant for lower scales because of the joint effect of the Kolmogorov power 
spectrum that suppress contributions on scales s as | 33 | and the mixing efficiency that is suppressed as 

A Kolmogorov power spectrum on scales between 1 kpc and 10 kpc is then used in the following for the simulations 


of the turbulent magnetic field with the method described in Sec III 


B. Chandra observations of Hydra A 

Observations of the Hydra A galaxy cluster with Chandra |3S] are analyzed following the method outlined in HZ]. 
The source has been observed in 1999 and 2003 IMlISZj with the AGIS instrument composed of two arrays of imaging 
GGDs that are sensitive to X-rays between 0.2 and 10 keV |3S|. The energy resolution of AGIS of 0.1 keV at 1.5 
keV |39| enables the search for spectral irregularities due to photon/scalar mixing. A total live-time of 238 ks is used 
in the analysis. 

The spectral analysis, performed with the XSPEC package version 12.7.1, shows that the spectrum of the cen¬ 
tral source is well-modeled (x^/rxdof = 48.59/58) by a power-law with a photoelectric absorption term [dN/dE = 
e~'^) avec F = 1.52 ± 0.17. The optical depth r is due to photoelectric absorption and is described by the 
hydrogen column A^h = 2.54 ± 0.33 x 10^^ cm“^ as r = iVH(T(if(l -|- z)) where cr{E) is the photoelectric cross-section 
as a function of the energy. As the absorption is mostly due to the intracluster medium of Hydra, the redshift of the 
source z = 0.0538 is taken into account. The spectrum is shown on Fig. [^together with the model and the residuals 
of the fit. Because of the strong photoelectric absorption, the central source is not visible above 1 keV and the search 
for spectral irregularities is therefore restricted to energies higher than this threshold. 


C. Statistical analysis of the irregularity level in the X-ray spectrum 

In order to estimate if some spectral irregularities induced by the mixing between photons and the scalar field can 
be accommodated by the data on Hydra A, they are fitted to the spectrum obtained in the previous section. Assuming 
that the intrinsic spectrum of the central source is well modeled by a power law with an absorption term (see dN/dE 
above), the data are fitted with this spectral shape modulated by the irregularities due to the scalar-photon coupling. 
The spectral shape fitted is then dN/dE x where is obtained by the propagation of an unpolarized photon 
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Energy [ keV ] 


FIG. 6: X-ray spectrum of the central source in Hydra A galaxy cluster observed by Chandra. Top panel: reconstructed 
spectrum with model (see text for details). Bottom panel: residuals to the model. 


beam in a random realization of the magnetic field. When A —>■ oo, P-y^.y = 1 and the initial spectral shape dN/dE 
is retrieved. The assumption of this particular spectral shape does not bias the result because a power law is the 
simplest spectral shape (in this context, the most regular) that could be used to describe the data. Assuming a more 
complicated spectral shape could yield a better constraint as the model would provide additional degrees of freedom 
and potentially accommodate artificially statistical fluctuations, thus excluding spectral irregularities more efficiently. 
As in ini. the goodness of fit is computed with a likelihood estimator C that assumes that the number of events in 
each reconstructed energy bin i follows a Poisson distribution: 

Si\ Bi\ ’ ^ ^ 

where in each reconstructed energy bin i, Si is the number of measured events from the signal region defined by a 
circle of V around the central source and Bi is the number of measured events from the background region defined 
as an annulus between 1" and 2.5" around the same position. The exposure rates of the two regions, corrected from 
their surface ratio, is written tg and tf, for respectively the signal and the background region. The expected signal 
photon rate rrii is predicted from the assumed spectral shape, convolved by the instrumental response functions: 

f AN ... 

m,= / TIF ^ P^^^{E)A{E)P{E\E)AEAE , (40) 

Jhin i Jo 

where A{E) is the true energy effective area and P{E\E) the migration matrix from true energies to reconstructed 
energies. The expected background photon rate is estimated from the measurement in the background region by 
requiring 


d\nC 

dbi 


= 0 , 


(41) 


giving the most likely value for bi. 

For each set of parameters {b, A), many simulations (1000) of the propagation of the system in the Hydra A magnetic 
field are performed. Each simulation corresponds to different realizations of the turbulent magnetic field and therefore 
gives different spectral irregularities patterns. To account for the lack of knowledge on the exact realization of the 
turbulent magnetic field, a likelihood ratio test is computed where the configuration of the magnetic field is considered 
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as a nuisance parameter |40| : 


X{b,A) 


sup£(6, A, 0) 
e 


sup C{b, A, 9) 
b,A,e 


(42) 


The log-likelihood ratio test —2 In A follows a distribution with two degrees of freedom |40] and the exclusion 
obtained on (5, A) at 95% C.L. corresponds to the set of parameters {b,A) for which —21nA(5, A) > 6.17. The 
obtained constraints on the model parameters are presented in Sec. |vg 


V. UV SPECTROSCOPY OF PKS 2155-304 

Because of the expression of E^, it is possible to probe higher values of b than what can be achieved with X-ray 
observations by searching for spectral irregularities at lower energies. This is rendered possible by looking at lower 
values of A, which corresponds to lower values of E^. The spectroscopy of a bright optical-UV point-like souce 
embedded in a galaxy cluster is used to obtain the constraints. In this study, we use the spectroscopy of PKS 2155- 
304, which is a BL Lac object around which a galaxy cluster is observed. This object has already been used in m 
to search for ALPs via spectral irregularities at TeV energies. Here, we take over the description of the magnetic field 
in the galaxy cluster of PKS 2155-304 that was used in m, where the magnetic field was conservatively described 
with a strength of 1 fiG and a Kolmogorov turbulence power spectrum on spatial scales between 1 kpc and 10 kpc. 
A typical electron density in galaxy clusters of 0.01 cm“^ is assumed. 

The optical and UV spectra of PKS 2155-304 have been measured by the HRS instrument on-board Hubble. 
The ECH-B echelon grating was chosen for its high resolving power of 1:80000 that is required to identify spectral 
irregularities occurring on the smallest energy scales. Observations were carried out in May 1993, before the first 
servicing mission that installed the COSTAR optics. The spectrum measured between 259.4 nm and 260.6 nm with 
the exclusion of a feature between 259.97 and 260.07 is shown in figure]^ The spectrum is shown normalized to its 
average value of 7.4 x 10“^^ergcm“^s“^A . The fit of a constant to the data gives a of 718.1 for 416 degrees of 

freedom, indicating that the spectrum is affected by some spectral features. The method used to search for spectral 
irregularities induced by the coupling between photons and the scalar field is similar to the one used in the Chandra 
analysis. Here, an underlying constant spectral shape is assumed, and a x^ estimator is used to quantify the goodness 
of fit. While the presence of some spectral features in the measured spectrum may decrease the sensitivity to spectral 
irregularities, this method is conservative when deriving constraints. 


VI. CONSTRAINTS ON THE MODEL PARAMETERS 

The above-described observations are used to set contraints on the different parameters of the model. For the X-ray 
spectrum, the method described in section [IV C| based on Eq. is used to determine the confidence intervals of the 
fit and set constraints at the 95% confidence level (C.L.). For the UV analysis, as stated in Sec. |v] a analysis 
is performed. All constraints shown in this section are valid for arbitrarily low masses below 4 x 10“^^ eV. In this 
case, the critical energy E^ does not depend on the mass of the field and the constraint is therefore independent of 
this mass. The exclusion regions can be drawn in a A — M plane to display simultaneously the constraints on the 
conformal coupling scale and the disformal coupling scale. This is done in Fig.j^where the left panel shows a focus on 
low values of M below 5 x 10“® GeV and the right panel focuses on higher values. In these plots, the grey regions are 
excluded by the present analysis, and are valid for a matter density of 10“^^ g/cm^. As mentioned before, contraints 
from the star burning rate allow only M > 346 MeV. However these values being obtained for matter densities of 
the order of 100 g/cm^, we derive constraints for lower M values, which might still be relevant in the case of a lower 
matter density. The region between A = 10® GeV and A = 10^° GeV cannot be probed as it corresponds to the 
gap between X-ray and UV observations. When A is set to lower values, only the value of E^ lowers so that spectral 
irregularities are still present. These orders of magnitudes for the exclusions were anticipated in Fig. [^and Fig. 

Figure shows the obtained constraints in the A vs. b plane. Again, the grey areas are excluded, and are valid 
for masse m < 4 x 10“^^ eV and densities of 10“^"^ g/cm®. In this case, the sensitivity in A at low b is limited 
by the value taken by E^ and not by the statistical uncertainties on the spectral points, as the conversion always 
occurs in the efficient mixing regime in which the oscillation length is smaller than the magnetic domain size. When 
A goes to higher values, E^ is shifted to higher energies and no conversion occurs in the energy range of interest. For 
b ~ 10“®^, E^ and E^ are of the same order of magnitude. The sensitivity in A decreases roughly as b^ so that E^ 
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FIG. 7: Spectrum of PKS 2155-304 measured by HRS with the ECH-B echelon grating. 
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FIG. 8: Gonstraint on A and M at the 95% C.L. for two ranges of A values. The grey area shows the constraint which is valid 
for all masses lower than 4 x 10“^^ eV. 


is kept constant, otherwise it is shifted to lower values for higher b and then no conversion occurs. For b < 10“^^, 
^ El and the ALP case is retrieved and the corresponding constraints on A hold. 

We have summarised all these constraints in Fig. 10 where we have chosen the disformal scale M = y/mp\mG, 
which corresponds to models of massive gravity with a graviton of mass m^j. We have represented the results in the 
(A“^, Tog) plane mimicking the corresponding diagram for ALP’s where usually the coupling is the displayed quantity. 

First of all, very low masses m < Hq are the only ones allowed when A" ^ ^ 10 ^ GeV The UV exclusion zone 
that we have obtained in this paper overlaps with the laboratory results from PVLAS and light-shining-through-a-wall 
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FIG. 9: Constraint on A and b at the 95% C.L. The grey area shows the constraint which is valid for all masses lower than 
4 X 10“^^ eV. 


Massive gravity (M^ = mo mpi) 



Exciuded by optics experiments (10“” g/cm^) 


This work g/cm^) 


FIG. 10: Constraint on 1/A and m at the 95% C.L. for the modified gravity case. Exclusions obtained from the present analysis 
are shown together with those obtained in different matter density ranges. 


experiments (labelled optical experiments in Fig. 101. For lower values of A“^, our X-ray analysis excludes a zone in 
masses which overlaps with the possible masses obtained from the star burning rate and the disformal coupling of 
scalars to photons, and where the effect of the direct coupling of scalars to photons implies that A“^ < 10“^° GeV“^. 
The X-ray constraint is slightly more stringent than the star burning rate bound on A“^. All our new results from 
UV and X-ray spectra are only valid for masses below the plasma frequency ojpi ^ 4 x 10“^^ eV. We have also recalled 
on the same figure the fact that all these bounds have been obtained in different environments and they could be 
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relaxed in models where the disformal coupling becomes density dependent. 


VII. CONCLUSION 

In this study we have considered astrophysical effects of the coupling to radiation of low-mass scalar fields, in a 
way somewhat comparable to the study of axion-like particle interactions with photons. We have opened up the 
parameters space to include a disformal coupling, as it would appear in massive gravity models or scalar field-based 
alternatives to dark energy. Considering photons travelling through turbulent magnetic fields in galaxy clusters, we 
showed that these interactions could lead to chaotic absorption patterns in the spectra. Such features have been 
searched for in both UV and X ray bands, with null results. This allowed us to obtain constraints on the coupling 
scales A and M of the scalar to radiation, i.e. conformal and disformal couplings respectively. For A our constraints 
are complementary to laboratory experiments as we exclude regions between 10^ GeV and 10^^ GeV, where dedicated 
setups are sensitive to scales up to 10^ GeV. For M, our exclusion limits extend to arbitrarily high energies for certain 
values of A. They complement both LHC constraints (those are limited to couplings to baryons) and constraints 
from star burning rates (obtained in a very different range for matter densities). Our limits concern scalar fields with 
masses below 4 x 10“^^ eV. If M is interpreted in terms of massive gravity, our analysis allows testing graviton masses 
down to TO > 10“^^ GeV. This means we probe modifications of gravity on scales smaller than one parsec, which are 
not allowed as long as A < 10^^ GeV. These results partly confirm the exclusions obtained from the study of star 
burning rates in a totally different matter density range. For larger values of A, X-ray observations are not sufficient 
and one may envisage gamma-ray probes in order to further constrain the parameter space of photons coupled to 
light scalar fields. 
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